SEPARATRIX MAP ANALYSIS FOR FRACTAL SCATTERINGS IN 
WEAK INTERACTIONS OF SOLITARY WAVES 
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Abstract. Previous studies have shown that fractal scatterings in weak interactions of solitary 
waves in the generalized nonlinear Schrodinger equations are described by a universal second-order 
separatrix map. In this paper, this separatrix map is analyzed in detail, and hence a complete 
characterization of fractal scatterings in these weak interactions is obtained. In particular, scaling 
laws of these fractals are derived analytically for different initial conditions, and these laws are 
confirmed by direct numerical simulations. In addition, an analytical criterion for the occurrence of 
fractal scatterings is given explicitly. 
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1. Introduction. Solitary wave interactions are important phenomena in sci- 
ence and engineering [U [5] . These interactions can be divided roughly into two types 
depending on the strength of the interactions. Strong interactions, often called col- 
lisions, are the interactions of solitary waves at close distance. They would occur 
when two solitary waves are initially far apart but move toward each other at mod- 
erate or large speeds. Weak interactions are the interactions of solitary waves at far 
distance through weak tail overlap. These interactions would occur if the two waves 
are initially well separated, and their relative velocities are small or zero. In inte- 
grable systems, collisions of solitary waves are elastic [T], and their weak interactions 
exhibit interesting but still simple behaviors [H [3l IH O [6] . For certain integrable sys- 
tems perturbed by higher-order corrections, if they can be asymptotically transformed 
to integrable equations, then their solitary wave interactions would closely resemble 
those in integrable systems [71IH1IS]. If the systems are non-integrable, however, soli- 
tary wave interactions can be extremely complicated, and they can depend on initial 
conditions in a sensitive, fractal manner. This fractal-scattering phenomenon was 
first discovered for kink and antikink collisions in the 0'* model \10[ [TlJ [12] , later in 
several other physical systems as well [131 [T3]. For these strong interactions, a reso- 
nant energy exchange mechanism between the collision and internal/radiation modes 
was found responsible for this fractal scattering. To analyze these fractal scatter- 
ings, approximate collective-coordinate ODE models based on variational methods 
[15] have been derived, and these ODEs are found to exhibit qualitatively similar 
fractal scatterings as in the PDEs [TH [ini [13 [TH] . Goodman and Haberman further 
studied these collective-coordinate ODE models using dynamical systems methods 
[TOl [ini [m [2H HH]- Performing asymptotic analysis along separatrix (homoclinic) 
orbits, they derived separatrix maps which led to the prediction of n-bounce reso- 
nance windows. It is noted that the separatrix maps derived in [TH [20l HU [22l [28] 
contain parameters which depend on initial conditions. In addition, these maps differ 
from one PDE system to another. On weak interactions, fractal scatterings have been 
found as well in a weakly discrete sine-Gordon equation and a class of generalized 
nonlinear Schrodinger (NLS) equations [231 [531 [5S] . For these weak interactions, the 
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mathematical analysis can be made more rigorous and quantitative. Indeed, by ex- 
tending the Karpman-Solovev perturbation method [3] , Zhu and Yang derived a sim- 
ple and asymptotically accurate ODE model for weak interactions in the generalized 
NLS equations with arbitrary nonlinearities |25| . After various normalizations, this 
ODE system contains only a single constant parameter which corresponds to different 
nonlinearities in the PDEs. These ODEs are a two-degrees-of-freedom Hamiltonian 
system with highly coupled potentials, and their forms are quite different from the 
collective-coordinate ODE models derived and studied previously for strong inter- 
actions [ni[Il[T71[THl[Tni[lDl[lIl[12]. Zhu, Haberman and Yang further analyzed 
this ODE model and derived a simple second-order map by using asymptotic meth- 
ods near separatrix orbits [521 123 • ^ remarkable feature of this map is that it does 
not contain any free parameters after various rescalings, thus it is universal for all 
weak interactions of solitary waves in the generalized NLS equations with arbitrary 
nonlinearities. Despite its simplicity, this map can capture all the fractal-scattering 
phenomena of the original PDEs and the reduced ODEs very well both qualitatively 
and quantitatively [26j . Reduction of weak-interaction dynamics from the PDEs into 
a simple and universal second-order map is the main contribution of [25l [26l [27] . With 
the availability of this universal map, one may now expect a complete characteriza- 
tion and understanding of fractal scatterings in weak wave interactions. For instance, 
we now would like to know under what conditions fractal scatterings would occur or 
would not occur. For another instance, we now would like to know how these frac- 
tals change as nonlinearities of the PDEs and initial conditions of the solitary waves 
vary. In addition, we now would like to understand how the zoomed-in structures of 
the fractal are related to the original structures, and how these geometric structures 
dictate the features of interaction dynamics. All these questions can be answered by 
a careful analysis of this universal map. 

In this paper, we analyze this universal map in detail, which will provide a com- 
plete characterization of fractal scatterings in weak wave interactions in the PDEs. 
First we will show that this map has a fractal structure of its own. We will delineate 
the map's fractal by tracking its singular curves. Then we will connect the map's 
fractal to that of the PDE, and thus reach a deep understanding of the PDE's fractal 
as well as its solution dynamics. In addition, we will determine how the PDE's fractal 
changes when the soliton parameters and the nonlinearity of the PDE vary. A precise 
analytical criterion for the occurrence of fractal scatterings in the PDEs will also be 
given explicitly. All our analytical results are confirmed by direct numerical simula- 
tions. These results significantly advance our understanding of fractal scatterings in 
weak interactions of solitary waves. 

2. Previous work. First, we summarize previous relevant work which will form 
the basis for our later analysis. We consider weak interactions in the generalized NLS 
equations 

(2.1) iUt + U^^ + N{\U\^)U = 0. 
These equations admit solitary waves of the form 

(2.2) [/ = $(a; - <e)e*"*, 

where ^{9) is a localized positive function, ^ — Vt + xa is the wave's center position, 
(j) = V{x — £^)/2 + {(3 + V'^/4:)t — (To is the phase function, and /3 is the propagation 
constant (which determines the amplitude of the wave). This wave has four free 
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parameters: velocity V, amplitude parameter /?, initial position xq, and initial phase 
constant (Tq. In weak interactions, two such solitary waves are initially well separated 
with small relative velocities and amplitude differences. Then they would interfere 
with each other through tail overlapping. When time goes to infinity, they either 
separate from each other at constant velocities, or form a bound state. The exit 
velocity, defined as AVoo = IV2 — Vi |i^oo, depends on the initial conditions of the two 
waves. When the two waves form a bound state, we define A14o = 0. Throughout 
this paper, we label the left and right waves by numbers 1 and 2 respectively. 

In [5S], we have shown that for a large class of nonlinearities N{\U\'^), this weak 
interaction depends on the initial conditions in a sensitive, fractal manner. To ana- 
lyze this fractal-scattering phenomenon, an extended Karpman-Solov'ev perturbation 
method was utilized, and the following simple set of dynamical equations for soliton 
parameters were derived [25] : 

/2 3-) Crr = COsV'e'^, 

ipTT — + £) sinipe'' . 

Here 

(2.4) V' = A0, C = -V^Ae, r=^^^^^t, 

and A(j) are the distance and phase difference between the two waves, /3 — {(3ifi + 
/?2,o)/2, /3fc,o [k — 1,2) are the initial propagation constants of the two waves, c is 
the tail coefficient of the solitary wave with propagation constant /?, and P{f3) is the 
power function of the wave. This ODE system is universal for the PDE (|2.ip . and 
different nonlinearities N{\U\'^) only correspond to different constant parameter e. If 
e = [such as when Eq. (|2.1[) is the original NLS equation], the ODEs (|2.3p are 
integrable, and their solutions have explicit functional expressions [HKl?]. If e 7^ 0, 
(j2.3p is not integrable. In this case, when e > and under certain initial conditions, 
we have found that these ODEs exhibit fractal scattering structures which agree with 
those in the PDEs (|2.ip both qualitatively and quantitatively. But when e < 0, no 
fractal scatterings arise in these ODEs and their corresponding PDEs under any initial 
conditions 

In order to further understand these fractal scatterings, we have analyzed the 
ODE system (|2.3p extensively in for |e| <C 1, using perturbation methods near 

separatrix orbits. These ODEs are a two-degree-of- freedom Hamiltonian system with 
the conserved Hamiltonian 

(2.5) i7(C,C>>)=i? + TJ7T^^', 

2(1 + e) 

where 

(2.6) E^^iC^ -iP^)-e'^cosiP 

is called the energy. We also define the momentum AI of Eqs. (|2.3p as 

(2.7) M ^ (ip- e'^ sinV'. 

Both E and M are conserved when e = 0, but vary over time when e 7^ 0. If the orbits 
are escape orbits where Coo = Clr^oo — ~oo (i.e. the two solitary waves eventually 
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separate from each other after weak interactions), then the exit velocities |^|oo can be 
calculated from Eqs. (|2.5p and (|2.7p as 



(2.8) \tU^^H+./WTMUil + e). 

In addition, il^oo, which determines the amplitudes of exiting solitary waves, can also 
be obtained as "000 = -^oo/Coo. Thus Mao is a key parameter for the prediction of 
weak- interaction outcomes. In order to calculate M^o, we notice that for weak wave 
interactions and when |e| ^ 1, the ODE solutions [such as C(''')] oscillate near a 
sequence of separatrix orbits of the unperturbed (e = 0) system before they escape to 
infinity. On these separatrix orbits, E — M — 0. If we consecutively enumerate the 
minimums of C(r) (where interactions are the weakest) and denote their energy and 
momentum values as i?„ and M„ (where n is the index of the C-minimum), then we 
can analytically calculate En and M„ successively by integrating along the separatrix 
orbits of the unperturbed system. This was done in [SJ , and we found that for 
|e| <C 1, En does not change, i.e. En = -Eo for all n > 1. But M„ does change. In the 
asymptotic limit of En ^ 1, il/„ ^ 1 and Mn/En <C 1, we found that the change of 
Mn is asymptotically governed by the following second-order separatrix map 



(2.9) Mn+i = Mn - sgn(g„ 

(2.10) Qn+l=Qn + 2Mn+l, 



with initial conditions Ma and Qo, where Qa = —SqMo, So = 5'|t=o, and 



7rRe(C) ' ^=V^^' ^ = -c"'^°*^(^)- 

Here the multi-valued functions and acoth(.) are chosen uniquely by requiring 
Im(-y^) > and Im(acoth(.)) e [0, tt) at the initial time. The variable Qn in the 
above map is an auxiliary variable which is related to the function S. 

When the ODEs (|2.3p are integrable (e = 0), 5' is a conserved quantity. In this 
case, the integrable solution Q{t) develops finite-time singularity if 5*0 = 2k, Re(Co) ^ 
0, fc = 0, ±1, ±2, or Re(Co) = 0, Im(Fo) = 0, see [27] for details. These finite-time 
singularities play important roles in the formation of fractal scatterings. Indeed, it 
was observed from numerical simulations in [25 that fractal structures for e 7^ 
bifurcate out from points where integrable solutions develop finite-time singularities. 
This fact will be proved in this paper by the analysis of the map. 

The map (|2.9p - (|2.10p can be normalized into a very simple form. Let 

(2.12) G=^^, m„ = G-i/3M„, g„ = G-i/3Q„, 

TT 

then this map becomes 

/o1q^ sgn(£g») 

(2.13) run+i = mn 5 , 

(2.14) qn+i = g„ + 2m„+i. 

This normalized map is second-order, and it does not contain any free parameters 
(except a sign of e). Thus it is universal for all weak two- wave interactions in the gen- 
eralized NLS equations (j2.ip with arbitrary nonlinearities. For positive and negative 



Separatrix map analysis for fractal scatterings 



5 



signs of £, we will show this map has completely different behaviors. This explains 
why fractal scatterings appear in ODEs (|2.3p only when e > 0, but not when e < 0. 
Details will be given in Sees. 2] and [SI 

To demonstrate the validity and accuracy of the above simple separatrix map 
for describing fractal scatterings in the original PDEs (|2.ip . here we compare the 
fractal scattering structures obtained from the PDEs (|2.1[1 . the ODEs (|2.3p . and the 
map p.l3|) - ()2.14p . In all our comparisons below and in later sections, we take the 
cubic-quintic nonlinearity 

(2.15) N{\Uf)^\Uf + S\U\^ 

with 6 = 0.0003 in the PDE (|2.ip (comparisons with other forms of nonlinearities are 
similar, see [25]). We take two types of initial conditions for the two solitary waves in 
the PDE (HH). One is that 

(2.16) /3oi=/3o2 = l, 

where the two waves initially have the same amplitudes. The other one is that 

(2.17) /3o4 = 1.0325, /3o,2 = 0.9675, 

where the two waves initially have unequal amplitudes. In both cases, the other 
parameters in the two initial solitary waves are the same as 

(2.18) xo^2 = -xo,i = 5, Vo,i = Vo^2 = 0, 0oa = 0. 

That is, the two waves initially have equal velocities and are separated by 10 spatial 
units. The initial phase of the first wave (^o,i can always be set as zero by phase 
invariance of the PDE l|2.ip . thus it does not constitute a restriction on the initial 
conditions. For both types of initial conditions, A(/)o(= ^0,2) is used as the control 
parameter. Corresponding to both types of initial conditions, we find from Eq. (|2.4p 
that e = 0.001 in the ODEs (j^ . 

For the first type of initial conditions (|2.16p . the corresponding initial conditions 
of the ODEs are 

(2.19) Co = -10, Co = ^0 = 0, 

and ipo{— (j)Q^2) is the control parameter. For the map (|2.13p - (|2.14p . the corresponding 
initial conditions can be readily found to be 

(2.20) qa = mo = -^TT^/^ tanipo ■ e~'^^^ 

since 6*0 = —1 in view of the expression (|2.1ip for S. Here we have assumed cos (f>Q > 
since we have shown before [27] that fractal scatterings arise only in this case. 

For the second type of initial conditions (|2.17p . the corresponding initial condi- 
tions of the ODEs are 

(2.21) Co = -10, Co = 0, ^0 = -0.01167. 
The corresponding initial conditions of the map are 



(2.22) 
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where 

(2-23) T(V.o) ^ ^ 

2cosi/'o + ?/'oe~ 

and -S'(^/'o) is defined in ([2lT|) . 

For each of the above two types of initial conditions, we have simulated the PDE 
(f2?T|) . the ODEs ([O]) . and the map ([2A3| - ([2A4|) for various values of the control 
parameter A(j)Q, and obtained the corresponding exit velocities AVoo [in the case of 
the ODEs and the map, the exit velocities have been obtained and properly rescaled 
through the formula (|2.8[) and relations (|2.4p . (|2.12[) ]. In all cases, we observed fractal 
scattering phenomena. For the first type of initial conditions (j2.16p (where the two 
waves initially have equal amplitudes), the fractal structures from the PDE, ODE 
and the map are displayed in Fig. 12.11 Here the fractals are symmetric in A0O: 
and they lie in a narrow interval near the A(j)Q = point (only a segment of the 
parameter region is shown for better visualization). We see that the three structures 
match very well both qualitatively and quantitatively. For the second type of initial 
conditions (|2.17p (where the two waves initially have un-equal amplitudes), the fractal 
structures from the PDE, ODE and the map are displayed in Fig l2.2l Here the fractals 
lie in a wide interval of A(/)o e [0,7r], and they do not possess any symmetry. Again 
these structures match very well. These comparisons show that the intricate fractal 
scatterings in the PDEs (|2.ip and the ODEs (|2.3p CAN be accurately predicted by 
the simple map (j2.13p - (j2.14p . These are the main results which have been obtained 
before in [HI 111 [23. 

In order to reach a complete understanding and characterization of these fractal 
structures and their solution dynamics in weak wave interactions, the separatrix map 
(j2.13[) - (|2.14[) needs to be carefully analyzed. This map exhibits a fractal structure 
of its own in its initial-condition space (when sgn(e) = 1) [26]. If this fractal of the 
map is well understood, then through connections of variables between the map and 
the PDEs/ODEs, the fractals in the PDEs/ODEs will be well characterized. In the 
following sections, we will analyze this map and use this information to delineate 
weak-interaction dynamics in the PDEs/ODEs. 

3. Analysis of the fractal in the separatrix map. The separatrix map 
(|2.13p - (|2.14p depends only on the sign of e, which can be 1 or —1. It turns out that 
this map has completely different dynamics for sgn(e) = 1 and —1. When sgn(£) = 1, 
this map is fractal-bearing, while when sgn(e) = —1, it is not [26]. This is consistent 
with the numerical observation of f2F that fractal scatterings occur in the ODEs (|2.3p 
only when e > but not when £ < 0. In this section, we analyze the fractal in this 
map when sgn(£) = 1. 

For convenience, we rewrite the map (I2.13p - (l2.14p as 

(3.1) ( ^"+1 ^" 

V m„+i J \ m„ 

where 

It can be seen that is differentiable when q ^ 0. The Jacobian of T is 



(3.3) = 



2 sg"('?) 1 
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Fig. 2.1. Comparisons of exit-velocity fractals from (a) PDE simulations, (b) ODE predictions, 
and (c) map predictions for the equal- amplitude initial conditions i2.16\l . The marked segment in 
(b) is amplified in Fig. I5.il where the solution dynamics in this segment is also shown. Labels 
■0^ (1), i/)"^ (1, 1), . . . are locations of singularity peaks (of infinite height) in (c), which correspond to 
intersections of the initial-value curve with the map's singular curves 7(1), 7(1, 1), . ■ ■ in Fig. \4. 1\ 
(see Sec. \4.l\ f0r details). 



0.16 




Fig. 2.2. Comparisons of exit-velocity fractals from (a) PDE simulations, (b) ODE predictions, 
and (c) map predictions for the unequal- amplitude initial conditions \2.17}) . For the segment of 
structures marked in (b), its movement with varying values of e is displayed in Fig. \4-4\ Labels 
V'f , ■ • ■ '^^^ locations of singularity peaks in (c) which correspond to intersections of the initial- 

value curve with the map's singular curves 70,7(1), . . . in Fig. \4-3\ (see Sec. \4.2\ for details). 
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The determinant of this Jacobian matrix is equal to one, thus the map T is area- 
preserving and orientation-preserving. Chaos if it occurs wiU be Hamiltonian chaos 

m- 

It is easy to see that T'^{—q, — m) = -~!F"{q, m), thus JF" is anti-symmetric with 
respect to the origin. It can also be seen that T does not have any (bounded) fixed 
points, but it has a lot of periodic orbits. For example, {(2~^/^, 2~^/^), (— 2~^/^, — 2~^/^)} 
is a period-2 orbit, and {(1, 1), (1, 0), (— I, — I), (— I, 0)} is a period-4 orbit. One can 
easily verify that these two periodic orbits are both unstable. All the other periodic 
orbits are unstable too as we will see later in this section. 

An important property of the map is that it is reversible in M^\r2, where M is 
the set of real numbers, and 

n = {((j,m) e : g = 2m}. 
Indeed, J-^^ has an explicit formula when q ^ 2ra: 

(") r+ iisr ) ■ 

In the later text, we will utilize two subsets of the (g, m) plane which we define here: 

n+ = {(q, m) e : g - 2m > 0}, 

and 

VL" ^ {{q,m) <eM? : q - 2m < 0}. 

3.1. Singular curves of the map and their identifications. The limit be- 
haviors lim„^oo ^"'{qoi tt^o) a-re important as they correspond to the outcomes of weak 
interactions. Indeed, from moo — lim„^oo "nin and the scaling relations (|2.12p . we can 
obtain the exit velocities |Coo| from Eq. (|2.8p . For almost all initial points {qo,mo), 
rrioo exists and is finite (while ^oo = ±oo). Thus orbits of almost all initial points 
eventually escape to goo = c» or — oo along the horizontal (constant-m) direction. But 
there are two special sets of initial points which are different. One set is the initial 
points which make qn = for some n > 0, i.e. 

(3.5) S = {{qo, mo) e M2 : g„ = 0, n = 0, 1, 2, ...}. 

For these points, iterations can not continue for (^n+i, m„+i) since = 0, thus we 
call them singular points, and S the singular set. On these points, we formally let 
m„_|_i = oo, and consequently moo = oo as well. For singular initial points, the exit 
velocities of weak interactions are infinite [see Eq. (|2.8|) ]. thus they correspond to 
peaks (of infinite height) in Fig. I2.ir c). The other set is the initial points (go,"^o) 
which are periodic or quasi-periodic. On these points, m„ oscillates forever, hence 
moo does not exist. These points correspond to spatially- localized and temporally 
oscillating bound states in the PDEs (|2.ip . In this case, we set moo — 0, which gives 
zero separation velocities |Coo| from Eq. (|2.8[) (recall that < for fractal scatterings 
[27]). This way, moo can be defined everywhere in the initial value plane (go,"^o)- 
Numerically we have computed |moo| over this plane (iterating 500 steps instead of 
infinite steps), and the result is plotted in Fig. 13.11 (here color levels correspond to 
|moo| values). This graph is a fractal, as can be easily verified by repeated zooms into 
it. 
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Fig. 3.1. Values of \ 

TTioo I the initittl-vdlue space (qO)^o) the Tnap \2. 1 3\) - \2. 14\ with 
sgnis) = 1. Colors represent value levels of |moo|. 



In Fig. 13.11 singular points form an infinite number of smooth, eurves. Eaclr of 
these singular curves lies in the middle of a red stripe (of varying thickness), and 
all these singular curves form the backbone of the map's fractal in Fig. 13.11 These 
singular curves are directly related to the singularity peaks (of infinite height) in the 
map's exit- velocity fractals in Figs. I2.ir c) and I2.2r c). This is because on a singular 
curve, moo = oo, thus the corresponding exit velocity is also infinite [see Eq. (|2.8p ]. 
In view of this, these singular curves correspond to the singularity peaks in the ODE's 
exit-velocity fractals and counterpart structures in the PDE's exit-velocity fractals. 
These singularity peaks (or their counterparts) in turn form the backbones of the 
map's, ODE's and PDE's exit- velocity fractals in Figs. 12.11 and 12.21 Thus if we can 
clearly characterize the singular curves of the map, then a good understanding of the 
fractals in the ODEs and PDEs will be reached. In the rest of this section, we focus 
on these singular curves. We will determine where they are located, how to identify 
them, what dynamics they represent, and what their asymptotics are at large or small 
values of go ■ 

About these singular curves, each one is characterized by a unique finite binary 
sequence a = sgn(go, 9i, ?«)■ Singular points on the same curve have the same 
binary sequence, while different singular curves have different sequences. Let us denote 
the singular curve with a binary sequence a = (ao,ai, ■■■,an) as 7(a). Then 



(3.6) 7(a) = {(go, "lo) € : qn+i = 0, sgn((7o, 9i, g™) = a}. 
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Fig. 3.2. Plots of singular curves 7(09), 7(00 1 <Jl)i 7(iOi ii > C'2), 7(i0i ^ii C'2i C3), and 
7(ao, ai, a2, as, 04) for all binary numbers o/ao,...,a4. The dashed lines in the right plane are 
7(1), 7(1, 1), 7(1, 1, 1), 7(1, 1, 1, 1) and 7(1, 1, 1, 1, 1) respectively from left to right. The accumula- 
tion curve 7(1, 1, 1, . . .) is also shown (as dotted lines). Lines in the left half plane are similar. The 
box region will be magnified in Fia \.3..3[ 

Here we allow a to be empty, in which case we get the simplest singular curve 

(3.7) 70 = {((70, mo) : qq = 0}, 

which is the vertical axis. When n = (where qi = 0), we see from the map p.2p 
that 

(3.8) 7(1) - {{qo, mo) : too - -<Zo/2 + l/ql qo > 0}, 

(3.9) 7(-l) = {(qo,mo):TOo = -qo/2-l/gg, go < 0}. 

These two curves are plotted in Fig. 13.21 They are located in the middle of the 
thickest red stripes in the right and left half planes of Fig. 13.11 

Now we determine the relations between singular curves with different binary 
sequences. From our definitions, we see that if (go, too) € 7(00, a), then T{qo,mQ) G 
7(a), where a is any finite binary sequence and uq — ±1. In addition, if ao = 1, i.e. 
qo > 0, then qi — 2toi > in view of Eq. (|2.14p . thus J-{qo,nio) E $1+. Similarly if 
ao — —1, then J-{qo, Too) G fi^. As a result, we have 



(3.10) 
(3.11) 



.F(7(l,a))=7(a)n^^+: 
.F(7(-l,a))=7(a)f|r!- 




Fig. 3.3. Amplification of the box region in Fig. \3. 2\ The dash-doted lines are 7(1, —1), 
7(1, —1, —1), 7(1, —1, —1, —1) and 7(1, —1, —1, —1, —1) respectively from right to left. Accumulation 
curves 7(1, —1, —1, —1, . . .) and 7(1, —1, 1, 1,1,.. .) are also shown. 

Or written differently, we have 

(3.12) 7(l,a)=.F-i(7(a)n^^+), 

(3.13) 7(-l,a)=.F-i(7(a)nf^")- 

Tliese relations tell us that each singular curve 7(a) has two pre- image singular curves 
7(1, a) and 7(— 1, a) under the map T . In particular, 7(1) and 7(— 1) are the two pre- 
image singular curves of 70. 

Thus from the basic singular curve 70 given in Eq. p.7p . two singular curves 7(1) 
and 7(— 1) in Eqs. (|3.8p - (|3.9p are obtained. We can successively construct 7(a) for 
any finite binary sequence a by repeatedly applying the inverse map J-~^ on 79. The 
first few of these singular curves are plotted in Fig. 13.21 It is easy to see from the 
map p.2p that for any binary sequence a, 7(a) and 7(— a) are anti-symmetric to each 
other about the origin, i.e. 

(3.14) 7(-a) = -7(a), 

thus we only consider curves on the right half plane with a = (1, ai, a„) below. 
From Fig. 13.21 we see that 7(1), 7(1,1), 7(1,1,1), ... form the primary cascading 
sequence from the left to the right. This sequence accumulates to the limit curve 
7(1, 1, 1, 1, ...) which is plotted as a dotted line in Fig. 13.21 and can be seen in Fig. 
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13.11 This primary sequence of singular curves defines the overaU geometry of the 
map's fractal in Fig. 13.11 On the left side of each primary curve, there is a secondary 
structure whose infinite curves are very close to each other and thus visually show as 
one "thick curve" in Fig. 13.21 Each secondary structure lies near its corresponding 
primary curve. To probe these secondary structures, we zoom into the box region 
of Fig. 13. 2( which contains a segment of the secondary structure for the primary 
curve 7(1). The result is shown in Fig. 13.31 From this zoomed- in graph, we see that 
among numerous curves in this secondary structure, there is a secondary sequence of 
curves 7(1,-1), 7(1,-1,-1), 7(1,-1,-1,-1),... which cascades to the left. This 
secondary curve sequence defines the overall geometry of this secondary structure. 
Near each curve in this secondary sequence (and on its right hand side), there is a 
higher-order structure which can be probed by repeated zooms. Every time one zooms 
into a higher-order structure, the cascading direction of its higher-order sequence is 
reversed, and the side of the higher-order sequence relative to its associated curve 
(either left or right) is also reversed. The binary sequences for these higher-order- 
sequence curves are {(a, 1), (a, 1, 1), (a, 1, 1, 1), ...} when the sequence cascades to the 
right and hes on the right side of 7(a), and {(a, —1), (a, —1, —1), (a, —1, —1, —1) ...} 
when the sequence cascades to the left and lies on the left side of 7(a). Here a is the 
binary sequence for the associated curve of this higher-order sequence. 

Based on the above pattern, we can identify any singular curve with an arbitrary 
binary sequence. For instance, to identify the curve with a binary sequence a = 
(1, 1,-1,-1,-1, 1, —1), we first go to the curve 7(1) (see Fig. 13.21) . find its associated 
primary sequence (which lies on its right hand side), and pick out the first member 
of that sequence (not counting 7(1) itself). The picked curve is then 7(1, 1), see Fig. 
13.21 Next, we go to the secondary curve sequence of 7(1, 1) (which lies on its left) 
and pick out the third member of that sequence, which is 7(1,1,-1,— 1,-1). Next 
we go to the higher-order sequence of 7(1, 1, —1, —1, —1) (which hes on its right) and 
pick out the first member of that sequence, which is 7(1,1,-1,-1,-1,1). Lastly 
we go to the higher-order sequence of 7(1,1,-1,-1,-1,1) (which hes on its left) 
and pick out the first member of that sequence, which will be the singular curve 
7(1, 1, —1, —1, —1, 1, —1) that we are looking for. 

We have noted earlier in this section that the map T has a lot of periodic orbits. 
Then an interesting question is where these periodic orbits are located in the map's 
fractal in Fig. 13. 1[ and how these orbits are related to the above singular curves. 
Clearly every periodic point can not lie on a singular curve 7(a) with a finite binary 
sequence a in view of the definition p.6p of 7(a). But a periodic point can be viewed 
as being on a singular curve with an infinite binary sequence a whose digits are the 
signs of q of the successive iteration points which repeat with the same period as 
the periodic point (here a singular curve with an infinite binary sequence can be 
defined as the limit of the singular curve with a finite binary sequence). Using this 
viewpoint, we can understand where periodic points should be located in the map's 
fractal. For instance, in the period-two orbit {(2^^/^, 2^^/^), (—2^^/^, —2^^/^)}, the 
first point (2^^/'^, 2~^/^) lies on the singular curve with the binary sequence a = 
(1, —1, 1, —1, . . .), while the second point (-2^^/"^, — 2~^/^) lies on the singular curve 
with the binary sequence a^— (-1,1,-1,1,...). For another instance, in the period- 
four orbit {(1, 1), (1, 0), (— 1, — 1), (— 1, 0)}, the point (1,1) lies on the singular curve 
with a binary sequence (1, 1, -1,-1,1,1,-1,-1,.. .), and the point (1, 0) lies on the 
singular curve with a binary sequence (1,-1,- 1,1,1,— 1,-1,1,...), etc. The singular 
curves with these infinite binary sequences can be identified by the same scheme as 
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we detailed above, thus we can ascertain where these periodic points are located in 
the map's fractal. Obviotisly these infinite binary sequences of periodic points are 
infinitely close to their finite truncations, and points on the singular curves with 
truncated finite binary sequences have very different trajectories from those of the 
periodic points. Thus the periodic points of the map T are all unstable. 

From the above singular-curve identification scheme, we see that all the singular 
curves on the right half of the (goi'Tio) plane lie between two special accumulation 
curves, 7(1,1,1,1,....) and 7(1, — 1, — 1, — 1, — 1, ...) (see Figs. 13.21 and [3.3p . So does 
the fractal of the map on the right half plane as well (see Fig. 13. ip . In the whole 
plane, the map's fractal lies between the two accumulation curves 7(1, 1,1,.. .) and 
7(-l,-l,-l,---)- 

Similar to the above tracking and identification of singular curves, we can also 
track how regions in the (g, m) plane move under the map T . This is helpful for 
us to see how the orbit of an initial point moves in the (g, m) plane. Let us denote 
the regions above 7(1, 1,1,.. .) and below 7(— 1, —1, —1, . . .) as I?0: and the region 
between 7(1, —1, —1, —1, . . .) and 7(— 1, 1, 1,1,-. .) as Then we find that 

(3.15) T-\v^)^v^\^v^. 

Other regions such as T^^iT)^ can be obtained with the help of singular curves 
whose pre-images under the inverse map have been detailed above (notice that 
the singular curve 70 lies in the middle of the region Vi ) . Details will not be pursued 
in this paper. 

3.2. Asymptotic behaviors of singular curves. Next we determine the asymp- 
totic behaviors of singular curves as go and ±cxd. These asymptotic behaviors 
will be needed for our derivation of scaling laws of the exit-velocity fractals in the 
ODEs and PDEs as e — > under the unequal-amplitude initial conditions ()2.17p . see 
Sec. 14.21 But they will not be needed for scaling laws of fractals with equal-amplitude 
initial conditions (|2.16p . see Sec. 14.11 Since 7(a) and 7(— a) are anti-symmetric to 
each other, we only consider curves on the right half plane with a — (1, ai, a„) 
below. 

The asymptotic behaviors of singular curves 7(1,01, ...,a„) as 0+ and +00 

can be obtained from 7(1) in Eq. (|3.8p and the recursion relations (|3.12p - (|3.13p . and 
we get the following results: 

1. For any singular curve 7(1, a), if (goj^^-o) G 7(li3-)j then 

(3.16) mo = g(72 + Ca + ..., 90 ^ 0+, 

where (2Ca, Ca) is the intersection point between 7(a) and 17. 

2. If (qot'tio) G 7(1j0.i, ...,a„) with = ... = a„ = 1 (i.e. on a primary singular 
curve), then 

(3.17) mo ^ -^^-J-^qo + DnQo"^ + qo ^ +00, 
where 

n + 1 

(3.18) D„ = l + Da = l. 

n 

3. If (go, mo) G 7(1, ai, a„, — 1, a) where ai = ... = a„ = 1 and a is an 
arbitrary finite binary sequence, i.e. when (qQ^mo) lies in the secondary 
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structure of a primary curve 7(1, ai, a„), then 

(3.19) +(!i±|^,-3/2^^^^^ 

where (2Ca, Ca) is the intersection point between 7(a) and fi. 

Notice from p.l7p and (|3.19p that for any binary subsequence a, 7(1, ai, .., a„, —1, a) 
tends to 7(1, ai, a„) as go ~* +00, i.e. the distance between them goes to zero 
at large go values. This explains why we could (and should) treat singular curves 
7(l,a, — l,a) with arbitrary binary subsequences a as the secondary structures asso- 
ciated with the primary curve 7(1, a) earlier in this section. 

To prove the first asymptotic result (|3.16p . let 



sgn(gi - 2mi) , 
(gi-2mi)2 ' 



(3.20) {qa,ma) = T ^(qi, toi) = (gi - 2mi, mi + 
Then using (|2.14p . we get 

«i go , sgn(go) 

(3.21) = + 

For {qo,mo) G 7(1, a), in view of the recursion relation (|3.12p . we see that (gi,mi) € 
7(a) nrj+. When qo 0+, ([3:201) shows that qi - 2mi 0+. Thus (gi,mi) ap- 
proaches the intersection point between 7(a) and Q, i.e. qi — > 2Ca and mi — > Ca. 
Substituting these asymptotic results into p.2ip . (|3.16p is then obtained. 

The second and third asymptotic relations (|3.17p and (I3.19P can be proved by 
the induction method and by utilizing the recursion relation (|3.12p . These proofs are 
elementary and will be omitted. 

4. Scaling properties of fractals in the ODEs and PDEs. We have known 
that fractal structures in the ODEs (|2.3p and PDEs (|2.ip are completely determined 
by the map's fractal, together with the initial-value connections (|2.4p and variable 
scalings (|2.12p between the map and the ODEs/PDEs. Utilizing the knowledge we 
have gained on the map's fractal in the previous section, we can now obtain a deep 
understanding on the fractals in the ODEs and the PDEs. This will be demonstrated 
in this section. For definiteness, we will use the two types of initial conditions (|2.16p 
and (|2.17p as examples. In both cases, as the initial phase difference tpo changes, the 
corresponding initial condition of the map forms a parameterized curve in the (go, "^o) 
plane. This curve, denoted as A"^, intersects the map's |moo| fractal in Fig. 13.11 and 
this intersection then completely determines the exit-velocity fractals of the ODEs 
and PDEs shown in Figs. OandO 

4.1. The case of equal-amplitude initial conditions. In this subsection, 
we consider the first type of initial conditions (|2.16p where the two solitary waves 
initially have equal amplitudes. In this case, the corresponding curve of the map's 
initial points in the (go, "^o) plane can be seen from (j2.20p as 

7^1/3 

(4.1) = {(go, mo) : go = mo = -^^jyr^ tani/'o, '4'o G (-7r/2, 7r/2)}. 




Fig. 4.1. Curve A*^ of the map's initial points 
initial conditions (2.1 (It . Some singular curves are 
those in Fig. \3. 2\ and the solid curves are 7(1, — 1) 



[4.1^ uihich corresponds to the equal- amplitude 
also shown. The dashed curves are the same as 
and 7(— 1, 1). 



The reason for the restriction ipo G (— 7r/2,7r/2) is that fractal scatterings can only 
arise in the ipo intervals where H < [27]. In the present case, H < corresponds 
to the interval ipo € (— 7r/2, 7r/2). This parameterized curve A*^ is a straight line with 
slope one in the {qo, mo) plane, see Fig. 14.11 

The map's |moo| fractal in the initial-condition plane (qQ^mo) (see the previous 
section) is very instrumental for the understanding of exit-velocity fractals in the 
PDEs and ODEs. First of all, from the map's |moo| fractal, together with the initial- 
value curve (|4.ip . the formula (|2.8p and various scalings, we can easily construct the 
map's exit-velocity fractal in Fig. I2.1f c). This exit-velocity fractal of the map can 
be readily understood. For instance, let us denote the ipQ value at the intersection of 
with a singular curve 7(a) as i/j^ia). At each ip^{a.), the map's exit- velocity graph 
has a singularity peak of infinite height. To illustrate, a few simple iP^{sl) values are 
marked in Fig. I2.ir c). Using this connection, the map's exit- velocity graph can be 
completely understood from the map's |moo| fractal. Then the exit- velocity fractals 
in the ODEs and PDEs can be similarly understood. To be specific, wc find that 
the primary window sequence in the exit-velocity fractals of Fig. 12. 1[ which cascades 
to the left with the first member being the widest window, are associated with the 
primary sequence of singular curves {7(1), 7(1, 1), . . .} (at the intersection with the 
set A"^). The dense secondary structure on the right hand side of each primary window 
in the exit- velocity fractals of Fig. 12.11 corresponds to the secondary structure of each 
primary-sequence curve in the map's |moo| fractal. In particular, the value ip^il, — 1), 
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which is marked in Fig. I2.ir c). corresponds to the secondary singular curve 7(1, —1) 
below the primary curve 7(1) in Fig. 13.31 If we zoom into each secondary structure 
of a primary window in the exit-velocity fractals of Fig. 12.11 we will see secondary 
window sequences which cascade to the right, i.e. the cascading direction of secondary 
window sequences is reversed from that of the primary window sequence. The reason 
for this is that in the map's \moo\ fractal (see Fig. 13. 2p . the cascading direction of 
secondary sequences of singular curves is reversed from that of the primary sequence 
as we have explained before. The exit-velocity fractals of the PDE, the ODE and the 
map in Fig. 12.11 can be zoomed further, and all their microscopic structures can be 
inferred from the map's singular curves in Fig. 13. 2i or from the map's |moo| fractal in 
Fig. I3.1l in general. One may notice that singularity peaks in the exit- velocity graphs 
appear only for the map and the ODEs, but not for the PDEs (see Fig. 12. ip . Near 
such singularity peaks, the two solitary waves collide and coalesce, which makes our 
reduced ODE model (|2.3p invalid. This explains the difference in those regions of the 
exit-velocity graphs between the ODEs and the PDEs. 

From the map, we can obtain another important piece of information on the 
PDE/ODE's fractal as the initial solitary-wave separation A^o (i-e- Co) varies. In the 
present equal-amplitude initial conditions (|2.16p and (|2.18p , if A^o takes other (large) 
values, we can easily see from (j4.1l) that the values of — rriQ are independent of 
A^o ■ This means that the exit- velocity fractal of the map [see Fig. I2.ir c)] will remain 
the same for different initial solitary-wave separations, which in turn implies the same 
for the exit-velocity fractals in the PDE/ODEs. This is a surprising fact, and it has 
been confirmed by our direct PDE/ODE simulations. It is noted that this fact does 
not hold for the unequal-amplitude initial conditions (|2.17p because {qo,mo) in such 
cases will depend on A^o, see (|2.22p . 

In addition to the above qualitative understanding of the exit-velocity fractals 
in the PDEs and ODEs, we can further obtain the scaling properties of these frac- 
tals, i.e. we can determine quantitatively how the fractal structures in the PDEs 
and ODEs change as the parameter e varies. For this purpose, we notice that the 
map's Imool fractal at the intersection with A*^ on the right half plane lies between 
two accumulation points, {qa,qa) = (0.741,0.741) and {qb,qb) — (1.271,1.271) on 
7(1, — 1, — 1, — 1, ...) and 7(1,1,1,1,...) respectively (see Fig. 14. ip . In view of the 
initial-condition connection (|4.ip . the corresponding -00 values of these two accumu- 
lation points are 

(4.2) Vl = -atan(2(767r-i/3£i/3)^ ^ ~atan(2ga7r-i/3ei/3). 

These ipL and V'fl, values are the left and right boundaries of the map's exit-velocity 
fractal on the negative A(/)o axis [see Fig. I2.1f c)]. and they are the map's predictions 
for the fractal regions in the PDEs/ODEs. These formulae show that V'L.i? as e —s- 
0+, which means that the whole fractal region shrinks to V'o = as e ^ 0+. Notice 
that when e = 0, the ODE solution C(r) under the equal-amplitude initial conditions 
(|2.19p develops finite-time singularity at V'o = 0, where Re(Co) = 0,lm(i^o) = (see 
Sec. [2]). Thus when e 0+, the fractal region shrinks to the ipo point which develops 
finite-time singularity in the integrable ODEs (see also [25]). Formulae (|4.2p further 
show that this shrinking is at the rate of e^/'^. To confirm this analytical prediction, we 
directly computed the exit-velocity fractals of the ODEs under the equal-amplitude 
initial conditions ()2.19p as e takes on smaller and smaller values of 0.1, 0.01, 0.001 and 
0.0001, and the results are displayed in Fig. 14.21 We see that as e — > 0+, the ODE's 
fractal region indeed approaches ipo = [see Fig. 14.2( 1-4)]. In addition, the fractal 
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Fig. 4.2. Left: change of fractal structures in the ODEs i2.3V as e ^ 0+ under the equal- 
amplitude initial condition fOtl) . The e values are (1) 0.1; (2) 0.01; (3) 0.001 and (4) 0.0001 
respectively. Labels ipL and ipn represent the left and right boundaries of the fractal region. Right: 
tan(i/'i) (blue) and ta.n{ipji) (red) versus e^^^. The circles are data from the ODE simulations, and 
the straight lines are the analytical formulae from the map. 

region's left and right boundaries tpL and V'i? indeed shrink in proportion to e^/^. 
Furthermore, the constants of proportion match the analytical values in Eq. (|4.2p as 
well [see Fig. 14.2( 5)]. Similar agreement has also been found for PDE fractals, see 

m- 

The comparison in Fig. 14.2( 5) indicates that for the present equal-amplitude 
initial conditions, the map's predictions are asymptotically accurate as e ^ 0+. This 
is not surprising, as the fractal region here lies near -00 = when |e| <C 1. In this 
case, we can easily see from the definitions (|2.6p and (|2.7p that <^ 1, Mq <C 1 
and Mq/Eq ^ 1 for the present initial conditions when ipo lies inside the fractal 
region. Thus the assumptions for the derivation of the map (|2.9|) - (|2.10[) are satisfied, 
and consequently the map's predictions are asymptotically accurate. For unequal- 
amplitude initial conditions (I2.21|) . however, the assumption Mq/Eq ^ 1 will not be 
met in general, thus the map's predictions will not be asymptotically accurate (even 
though they are still qualitatively accurate), see the next section for details. 

4.2. The case of unequal-amplitude initial conditions. Now we consider 
the second case of unequal-amplitude initial conditions ()2.17p . In this case, the cor- 
responding curve of the map's initial values in the (go,'7io) plane can be seen from 
([2:^ as 

(4.3) = {(go, mo) : <zo = -r(0o)5(Vo)e~'/', mo = T(Vo)e"'/', ^0 e [0,27r]}, 

where T{ipo) is defined in (|2.23p . S{ipo) given by (|2.1ip . and the other involved pa- 
rameters specified in (j2.21|) . Here H < in the entire interval of V'o G [0, 2tt], thus no 
restriction on i/jq is needed [27]. These curves at two e values £1 — 0.01 and £2 — 0.005 
are displayed in Fig. 14.31 Each is a closed curve, and it intersects a singular curve 
7(a) twice. Let us denote the ipo values at the two intersections as ipHs.) and ^'2(^)5 
with ipi{a) < "02 (^)- These ipf 2(3-) values are the singularity peaks of infinite height 
in the exit-velocity fractal of the map [see Fig. I2.2f cl]. On the singular curve 70, the 
intersection points would be denoted as -0f and 1^2- As e ^ 0+, these intersections 
are such that qq +00 in the right half plane and qq 0^ in the left half plane. 
This behavior can be readily understood by examining the intersections of A"^ with 
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the horizontal and vertical axes. The intersections with the horizontal axis are such 
that mo = 0. In view of (14. 3p as well as the expression (|2.23p for T{iPq), we see that 
these intersections occur when T{iPq) = 0, i.e. when t/jq — and tt, regardless of the 
s values. Note that under the present initial conditions, 

(4.4) qA = - lim r(^o)5(V'o) = 0.962, qb = - lim T(^o)5(V'o) = 0.641, 

thus these two intersections on the horizontal axis are {qA£~^^^ ,0) and {qB£~^^^ ,0), 
which move to (+cxd,0) as e — > 0"*". Similarly, the intersections of with the vertical 
axis occur when S'(V'o) = 0, i.e. when ^pA = 1.988 and -05 = 2.871, regardless of the e 
values. Then the two intersection points on the vertical axis are [0, r(-0yi)e^^/'^] and 
[0, T(?/'_B)e~"^/^], which approach (0, — oo) as e — > 0+. 

With the help of the map's singular curves as well as the initial-condition curve 
A*^ in Fig. 14.31 we can now understand the map's exit- velocity fractal in Fig. I2.2r c). 
and hence the PDE/ODEs' exit-velocity fractals in Fig. \2.2( a. b). For these initial 
conditions, e ~ 0.001 (see Sec. 2), and the initial- value curve A"^ goes outside the 
box of Fig. 14.31 (thus not displayed). Instead, we will use the curve A*^^ in Fig. 
14.31 (with £2 = 0.005) as a qualitative guide. The curve A^ in the first quadrant 
[above the accumulation curve 7(1, 1, . . .)] corresponds roughly to ipo G (tt, 2n). This 
segment of A^ does not intersect with any singular curves, thus its corresponding 
exit-velocity graph would be smooth (in Fig. 12.21 this segment of the graph is not 
shown). The right intersection point of A^ with 7(1, 1, . . .) is where ^^1(1: 1, . . .) ~ 0. 
This intersection corresponds to the left edge of the map's exit-velocity fractal in Fig. 
\2.2i c). From this intersection point down (leftward), A*^ passes through the primary 
sequence of singular curves (in the reverse order), which corresponds to the primary 
sequence of singularity peaks starting from A(j)Q « rightward (in the reverse order) 
in Fig. I2.2r c). The singularity peaks V'i(l) and "^Kli 1) in this primary sequence are 
marked in Fig. I2.2r c). which correspond to lower intersections of A^ with primary 
singular curves 7(1) and 7(1, 1). At the lower intersection of A^ with the vertical axis, 
■01 = V'A = 1.988, which is labeled in Fig. I2.2r c). From this lower vertical intersection 
point leftward, A^ passes through a thick band of singular curves, which corresponds 
to the structures right after the peak of 0f in Fig. I2.2r c) (these structures are not well 
resolved and only a few vertical points are visible). After this thick band, A"^ turns 
around (upward) and passes through another thick band of singular curves in the left 
half plane, which corresponds to the structures between the label '(c)' and the first 
major peak to its right in Fig. I2.2r c) (again these structures are not well resolved). 
The curve A^ crosses the vertical axis again (upper intersection) at -02 = V's = 2.871, 
which is labeled in Fig. I2.2f c). From this upper vertical intersection rightward, A*^ 
passes through the primary sequence of singular curves again (in forward order), 
which corresponds to the primary sequence of singularity peaks from V'f rightward 
and ending at 02(li 1, . . .) « tt in Fig. I2.2r c). From this correspondence between the 
initial- value curve A^ and the map's exit- velocity fractal in Fig. I2.2r c). a good and 
clear understanding of the exit- velocity fractals in the PDEs/ODEs [see Fig. I2.2r a. 
b)] is then reached. 

From the previous section, we know that when qq +00, the map's secondary 
singular curves below each primary singular curve approach this primary curve. In 
addition, it is easy to see that all singular curves in the left half plane approach 70 
(the vertical axis) when qq 0~ . Then in view of the small-e asymptotics of curve 
A^ described above, we see that as e ^ 0+, the intersections of A"^ with the map's 
singular curves approach the primary sequence {70, 7(1), 7(1, 1), 7(1, 1, 1), . . .}. As 




a result, when e 0+, the exit-velocity fractals in the PDEs/ODEs will converge 
to certain discrete -00 values whose corresponding {qo,mQ) points fall on the above 
primary sequence. As was explained in [27j . such discrete values are precisely the 
initial-condition points whose C-orbit in the integrable ODE system ()2.3|) develops 
finite-time singularities (this fact will be re-established again later in this section). 

In addition to the above qualitative descriptions of the PDE/ODE's fractals (see 
Fig. 12. 2[) . we can further determine how these fractal structures change quantitatively 
as £ ^ 0+. For instance, we can determine how the singularity peaks in the ODEs' 
exit-velocity fractal of Fig. I2.2r b) move as e varies [the PDE's fractal does not have 
singularity peaks but only counterpart structures, see Fig. I2.2f a)]. This can be 
done because we know how the curve moves with e [see Eq. (14. 3p and Fig. 14. 3j . In 
addition, its intersections with singular curves of the map have either qo — > -|-oo (in the 
right half plane) or qq — > 0~ (in the left half plane), where the asymptotics of singular 
curves have been derived in the previous section. Thus when the (go, mo) relation of 
the curve in (|4.3p is inserted into the asymptotic equations of singular curves and 
all variables are expressed in terms of the control parameter -00: quantitative changes 
of singularity peaks in the ODEs' exit-velocity fractals will be obtained. After simple 
algebra, we get the following asymptotic results in the limit of e 0+: 



(4.5) 



01 =0^, 1p2^1pB, 
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(4.6) 



.0(l,a)+Afc(a)e + ... 



(4.7) 



^1(1, a, -1, a) = V2(l, a) + Gfc(a)ei/2 + ^^(a, a)e'/^ + ... 



(4.8) 



M-h a) = V° + Pfe£i/2 + Qfe(a)e5/6 + . . . . 



Here a = (ai,...,a„), ai = ... = a„ = 1, a is an arbitrary finite binary sequence, 
Ak , Gfc , , Qk are constants whose values depend on the binary sequences behind 
them, Pk are a-independent constants, and fc = 1, 2. If n = 0, then a is empty (which 
is allowed). Eq. (|4.5p indicates that the singularity points i/jf, in the ODEs' fractal 
[see Fig. I2.2f b)] are e-independent. Actually these ipf, points in the ODEs/PDEs 
do depend weakly on e, but this weak e-dependence can not be captured by our 
map since it is beyond the asymptotic validity of the map. Eq. (14. 6p describes how 
a primary singularity point approaches its integrable counterpart as £ —^ 0"^, and 
this convergence is at the uniform rate of 0{e). Relations (|4.7|) and (|4.8p describe 
how the secondary structure of a primary singularity point approaches the integrable 
counterpart of this singularity point, and this convergence is at the uniform rate of 
0(e^/^). The length of the secondary structure, on the other hand, shrinks at the 
rate of 0(e^/®). By substituting the {qo,mo) relation of (|4.3p into the asymptotic 
equations p.lTp of primary curves and taking the limit of e ^ 0, we see that 



In addition, at these and ?/'^(l,a) values, it is easy to check that Re(Co) ^ 0, thus 
solutions of the integrable ODEs develop finite-time singularities in ({t) (see Sec. 2 
and [27j ) . Thus relations (|4.7p and (|4.8p , together with Eq. (|4.2p and the discussions 
below it, quantitatively prove that when e — )■ 0+, the fractal regions in the non- 
integrable ODEs shrink to the points which develop finite-time singularities in the 
integrable ODEs, as was originally observed numerically in [25]. 




Now we compare the above analytical e-scaling laws for the exit-velocity fractals 
in the PDEs/ODEs with direct numerical simulation results. Comparisons with only 
ODE simulations will be performed, as PDE simulations at very small e values are very 
expensive and time consuming. For the ease of comparison, we take a small segment 
of the exit- velocity fractal in the ODEs as marked in Fig. I2.2f b). which corresponds 
to the segment on the initial-value curve containing the lower intersection with 
7(1) and its secondary structure. Then we monitor how this segment of the fractal 
moves as e varies by directly simulating the ODEs (|2.3p . At three e values of 0.002, 
0.001 and 0.0005, these segments of the fractal structures in the ODEs are displayed 
in Fig. 14.41 (1-3) respectively. Here the vertical solid lines mark the singularity point 
in the integrable ODE system, the vertical dashed lines mark the singularity 
point V'^(l) in the non-integrable ODEs (with non-zero e), and labels ip1 ^(1) mark 
the left and right ends of the secondary structure. It is seen from these figures that 
as e — > 0"'', both and the secondary structure approach the singularity point 

tp'^il) of the integrable ODEs (as predicted). To determine the convergence rates, the 
graphs of — V'"(l)i i'l ni^) ~ '/'°(1) and i>%{X) — ^£(1) various values of e 

are plotted in Fig. I4.4r a. b, c) respectively. It is seen that as e ^ 0+, V''^(l) — ^'"(1) 
approaches a linear function in e, confirming the 0(e)-convergence of the primary 
singularity point ■0'^(1) toward the integrable singularity point '\lf'(X) [see formula 
(|4.6p ]. Quantities ^(1) — '0°(1) approach a linear function in e^/^, confirming the 



(4.9) 



5(^°(l,a)) = 2(n + l), 5(^D = 0, A: = 1,2. 
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Fig. 4.4. Left: segments of the ODE's fractal structures near the lower intersection between A'^ 
and the primary singular curve 7(1) under unegual- amplitude initial conditions I2.21\l [this segment 
for e = 0.001 was marked in the whole fractal structure of Fig. 1 2. Bf b )]. The e values are (1) 
0.002; (2) 0.001 and (3) 0.0005. Right: (a) - ^"(l) versus e; (b) V>£(1) - (blue) and 

ip%{l) - (red) versus e^l'^; (c) iJ%{l) - i>l{l) versus e^/e. 



0(e^/^)-convergence of the secondary structure toward the integrable singularity point 
[see formula (|4.7|) ]. In particular, it is seen that as e — > O"*", tp^i^) — ^/'"(l) and 
— V'^(l) approach linear functions in e^/^ with the same slope, confirming the 
formula (|4.7p that the coefficient of e^^^ is independent of the singular curves inside 
the secondary structure. The quantity ~ "01 (1) approaches a linear function in 

e^/^, confirming the 0(e^/^) shrinking rate of the length of the secondary structure 
[see (|4.7p ]. Thus our analytical predictions (|4.6p - (|4.8p on the order of convergence of 
ODE fractals as e — * 0+ are confirmed. 

Quantitatively, the coefficients in front of the order of convergence in the analyt- 
ical formulae (|4.6p - (|4.8|) do not match numerical values though. For instance, in Fig. 
I4.4f a). the slope of — with respect to e is found numerically to be —3.89, 

which differs from the analytical value of A{1) — —5.12 in formula (|4.6p . This means 
that in the present case of unequal-amplitude initial conditions ()2.2ip . our analytical 
formulae ()4.6p - ()4.8p are not asymptotically accurate when £ ^ 1, which contrasts the 
equal-amplitude initial-condition case in Sec. 14.11 (where our analytical predictions 
were asymptotically accurate). The reason for this is that our map (|2.9p - (|2.10p . or 
(|2.13p - (|2.14p . was derived asymptotically under the condition of M„/i?„ <C 1 p6l[27] . 
This condition was satisfied for the equal-amplitude initial conditions (|2.16p (see Sec. 
14. ip . but is not satisfied for the unequal-amplitude initial conditions (|2.2ip . Indeed, 
one can easily check that for ipo values in the fractal regions of unequal-amplitude 
initial conditions (see Figs. 12.21 and l4.4p . Mo/Eo does not tend to zero as e ^ 0, thus 
the map's predictions become asymptotically inaccurate as we have just observed. 
Qualitatively, the map's predictions are still correct as Fig. 14.41 has demonstrated. 

5. Dynamics of ODE and PDE solutions in the fractal structures. From 
the previous two sections, we have reached a clear and deep understanding on the 
fractal graphs of exit velocities in the ODEs and PDEs (see Figs. I2.1l and [2.2p . In this 
section, we explain the solution dynamics of the ODEs and PDEs on these fractals. 
Notice that these fractals consist of hills of various widths. The 'center' of each hill 
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corresponds to a singular curve in the map's fractal, thus each hill can be identified 
by the binary sequence of that singular curve. We will show that once the binary 
sequence of a hill is given, then the solution dynamics of the ODEs and PDEs on that 
hill can be ascertained. Throughout this section, sgn(e) = 1, where fractal scatterings 
occur. 

We first describe the ODE solution dynamics in the exit-velocity fractal. To 
demonstrate, we pick the hill of binary sequence a = (1, —1, 1) in the ODE's fractal 
for equal-amplitude initial conditions in Fig. l2.1f bV The segment which contains 
this hill is marked in that fractal. The amplification of this segment is shown in the 
inset of Fig. IS.ir a). where the widest hill is the one of this binary sequence a. To 
examine the dynamics of ODE solutions on this hill, we pick three points on the hill 
with ipa < i^b < ipc, where the middle point ipb is the singularity peak point, and the 
other two points are on the two sides of this singularity point. These three points 
are marked in the inset of Fig. IS.lf a). At these three points, the ODE solutions 
^(t) and ijj{t) are displayed in Fig. IS.lf a. b) respectively. Here the time t has been 
rescaled back to the physical time for easy comparison with the PDE dynamics below. 
We see that at the singularity peak point ipb, the C{t) solution oscillates three times, 
then develops finite-time singularity at ~ 408 and terminates there. Each local 
minimum of C is a saddle approach |27| where the separation between the two waves 
is the largest. Each local maximum of ^ is a 'bounce' point where the two waves are 
locally the closest and interact more strongly. The ip{t) solution is mostly flat, except 
that it exhibits spikes at the three bounce points whose sign sequence is (—1, 1, —1), 
which is opposite of the binary sequence a. This ipit) solution also terminates at the 
finite-time singularity time tc- At the other two points ipa and ipc on the two sides of 
the singularity peak ipb, the ODE solutions do not develop finite-time singularities. 
From Fig. IS.lf a. b), we see that these solutions are almost indistinguishable from the 
singular solution of ipb up to the singularity time tc- At t « tc, the ({t) solutions in 
both cases have a global maximum, where the two waves are the closest and interact 
most strongly. This time was called the collision time in |2Sj . Beyond this time, the 

solutions of ipa and ipc both go to — oo with finite speed. The ip solutions of ipa 
and Ipc approach constants as t ^ oo, but these asymptotic constants have opposite 
signs: the asymptotic constant for ipa at the left side of the singularity peak is positive, 
while that for ipc at the right side of the singularity peak is negative. At other points 
on this hill, the ODE solution dynamics is qualitatively similar to the ones above. 
From these examples, we can draw general conclusions for the ODE dynamics on a 
hill of an arbitrary binary sequence a = (ao, oi, 02, . . . , a„) in the exit-velocity fractal. 
At the singularity peak of the hill, the ((t) solution oscillates n -I- 1 times (i.e. the 
two waves bounce with each other n + 1 times), and then approaches infinity at a 
finite time tc and terminates. Here n -I- 1 is the length of the binary sequence a. The 
—ip solution exhibits spikes at the n + 1 bounce points before the singularity time tc, 
whose sign sequence is a. At other points of the hill, the ODE solutions are almost 
indistinguishable from this singular solution up to the singularity time tc- The time tc 
is approximately the collision time of all these ODE solutions. Beyond this collision 
time, all ({t) solutions go to —00, while all the ip{t) solutions approach constants 
which have the same sign on the same side of the hill but opposite sign between the 
two sides. 

From the above description of ODE dynamics on hills of exit-velocity fractals, we 
see that the 'physical' meaning of the binary sequence a of a hill in the ODE solutions 
is that a gives the sign sequence of the —ip{t) solution at the bounce points before the 
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collision time. In addition, the difference in solution dynamics between the left and 
right sides of the hill is that the asymptotic constants of their — solutions have 
opposite signs. These two facts can be readily explained. First, from the definition 
(j3.6p . we know that the binary sequence a is the signs of (go, 9i, ■ • ■ , 9n), with qn+i = 0. 
This orbit of the map corresponds to the ODE solution at the singularity peak of 
binary sequence a in the exit- velocity fractal. From Eqs. (|2.9p and p.l2p . we see that 
qk has the sign of —AMk, where AMk = Mk+i — Mk, and Mk is the momentum value 
at the fc-th saddle approach. Hence a is also the signs of (— AMq, — AMi, . . . , — AM„). 
Furthermore, from Ref. [57] [see Eqs. (3.4) and (5.6) in particular], sgn(— AM^) is 
the sign of —tp at the k-th ^-maximum (bounce point). Thus the binary sequence a 
is equal to the sign sequence of —-0 at the bounce points (before the singularity time 
or the collision time). Regarding the signs of — ^/'(oo) on the two sides of the hill, we 
notice from the map (|2.13p - (|2.14p that in the map's |moo| fractal in Figs. 13.11 and 
13.21 when {qo,mo) is on the left (right) side of the vertical axis 70 and in its vicinity, 
goo is positive (negative). Using the recursive relations p.l2p - p.l3p between singular 
curves and the orientation-preserving property of the inverse map we see that 

for {qo,mo) lying on the left (right) side of every singular curve (and in its vicinity), 
goo is positive (negative). Notice that like bounce points above, the sign of goo is the 
same as that of —■ip{oo). In addition, the two sides of a hill in the ODE's exit-velocity 
fractal correspond to the two sides of the singular curve in the map's |tooo| fractal. 
Thus the values of —ip{oo) on the two sides of a hill in the exit-velocity fractal have 
opposite signs. The specific signs of —tp{oo) on the two sides of the hill depend on 
how these two sides of the hill correspond to the two sides of the singular curve. If the 
left side of the hill corresponds to the left side of the singular curve, which is the case 
for hills starting from ■01 rightward in the unequal-initial-amplitude fractal of Fig. 
I2.2f b). then — 0(oo) would be positive (negative) on the left (right) side of the hill. 
But if the left side of the hill corresponds to the right side of the singular curve, which 
is the case for all hills in the equal-initial-amplitude fractal of Fig. I2.ir b) and the 
hills starting from ipf leftward in the unequal-initial-amplitude fractal of Fig. I2.2f b). 
then —ip{oo) would be negative (positive) on the left (right) side of the hill. 

The topic of higher interest to us is the interaction dynamics in the PDEs (|2.ip 
rather than in the ODEs (|2.3p . So next we describe the PDE solution dynamics in the 
exit-velocity fractal. This PDE dynamics can be predicted from the ODE dynamics 
above. To illustrate, we take two A0o values on the hill of binary sequence a = 
(1, —1, 1) in the PDE's equal-initial-amplitude fractal of Fig. I2.1f a). which correspond 
to the two 00 values of the ODEs on the two sides of the singularity peak as marked in 
the inset of Fig. IS.lT a). For these two A0o values, the PDE solutions are displayed in 
Fig. IS.lT c. d). We see that before the collision time tc ~ 408, the two PDE solutions 
are almost identical with each other. In this time period, the two waves bounce 
three times. At the three bounce points, the left wave has higher, lower and higher 
amplitudes sequentially (when compared to the right wave). At the collision time, 
the two waves are the closest and interact most strongly. Afterwards, they separate 
from each other and escape to infinity. The main difference between the two PDE 
solutions is that, at the left A0o value, the exiting right wave has higher amplitude, 
while at the right A0o value, it is just the opposite. These behaviors are common for 
most of the points on this hill in the PDE's fractal. The only exception is a small 
region in the middle of this hill which corresponds to the singularity peak point and 
its immediate vicinity in the ODE's fractal [see inset of Fig. IS.lf a)]. In that region, 
the C solutions in the ODEs are either infinite or very large at the singularity time 
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Fig. 5.1. Dynamics of the ODE and PDE solutions on the hill of binary sequence a = (1, —1, 1) 
in the exit-velocity fractals for equal- amplitude initial conditions I2.1ffi) . The inset in (a) is an 
amplification of the segment marked in the ODE fractal of Fig. \2.lY b}. where the widest hill is the 
one of binary sequence a = (1, —1, 1). Three locations ipa < tpb ^ ''Pc on this hill are marked by two 
triangles and a vertical dashed line, with the dashed line at ipi, being the singularity peak point. The 
ODE solutions at ipa, ipb ^.^d. ipc are plotted as dash-dotted lines (blue), dashed lines (red) and solid 
lines (black) in (a, b) respectively. The PDE solutions at points corresponding to ipa and ijic in the 
PDE fractal of Fig. \2.1V a) are shown in (c, d) respectively (contour plots). 



or collision time, which implies that the two waves collide and coalesce. When this 
happens, the reduced ODE model (|2.3p and its predictions become invalid. Indeed in 
the PDE's fractal, the central part of every hill dips down, which contrasts with the 
ode's fractal where the central part of each hill rises up to a peak of infinite height. 

The PDE solutions in Fig. IS.lT c. d) directly correspond to the ODE solutions at 
the two sides of the singularity peak in Fig. IS.lT a. b). In particular, the sign sequence 
a = (1, —1, 1) of the ODE's — solution at the three bounce points directly implies 
the "higher, lower, higher" amplitudes of the left wave at the three bounce points in 
the PDE solution, and the negative (positive) —ip{oo) values at the left (right) side 
of the singularity peak directly implies the lower (higher) amplitude of the exiting 
left wave in the PDE solution. These connections can be readily explained. Let us 
recall the relation A(j)t — A/3 as well as the fact that the propagation constant /3 is 
directly related to the wave amplitude [5S]. Then the sign of —ip, which is equal to 
the sign of — A/3 = /3i — P2, tells which of the two waves has higher amplitude. For the 



Separatrix map analysis for fractal scatterings 



25 



cubic-quintic nonlinearity (|2.15p . the amplitude is an increasing function of /3, thus 
positive —ijj means that the left wave has higher amplitude, which explains the above 
connections. 

Based on the above PDE examples and general ODE dynamics, we can draw 
general conclusions for the PDE dynamics on a hill of an arbitrary binary sequence 
a = (ao, ai, a2, . . . , an) in the exit- velocity fractal. For all points on the hill, the PDE 
solutions are almost identical to each other up to the collision times (whose values are 
almost the same for the entire hill). Before the collision time, the two waves bounce 
n + 1 times. At each bounce point, the left wave has higher (lower) amplitude if the 
corresponding digit in the binary sequence a is 1 (—1). Thus the physical meaning of 
the binary sequence of a hill in the exit- velocity fractal is that it gives the sequence of 
relative amplitudes between the two waves at the bounce points before the collision 
time, with digit 1 (—1) meaning the left wave is higher (lower). At the collision time, 
the two waves are the closest and interact most strongly. For most of the points on 
the hill (except a small section in the middle), the two waves separate from each other 
after the collision time, and the exiting waves have opposite relative amplitudes on the 
two sides of the hill. The specific relative amplitudes of the two waves are determined 
by the sign of —^'{oo) which has been given above. The left wave will have higher 
(lower) amplitude if ~'>p{oo) is positive (negative). 

6. The map for negative e and its predictions for the PDE/ODE sys- 
tem. In this section, we consider the e < case. For the cubic-quintic nonlinearity 
(|2.15p . negative e occurs when 5 < [25]. In this case, we will show below that fractal 
scatterings can not occur. The map p.l3p - (|2.14p for negative e differs from that for 
positive £ by only a sign, but this makes a crucial difference. Now the map is 



(6.1) J" 
and 



q — 2m 



(g-2m)2 



Notice that and J-'^^ are still area-preserving and orientation preserving, but the 
singular-curve structure changes drastically. To see that, we start with the singular 
curves 7(1) and 7(^1), which are now 

(6.3) 7(1) = {{qo,mo) : mo -qo/i - l/q^, qo > 0}, 

(6.4) 7(-l) = {{qo, mo) : mo = -go/2 + l/q^ qo < 0}. 



These two curves are shown in Fig. 16.11 We see that these curves are in the second 
and fourth quadrants. More importantly, they do not intersect with Q. Because of 
this, each of these two singular curves has only one pre-image curve under the map of 
J-'^^. This contrasts with the e > case where each singular curve has two pre-image 
curves. The pre-image curves of 7(1) and 7(— 1), 

(6.5) 7(1,1) =^-'(7(1)), 7(-l,-l) = ^-^(7(-l)), 



are also shown in Fig. 16.11 They do not intersect with Q either, thus have only one 
pre-image curve each as well. Repeating this process, then all the singular curves one 
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Fig. 6.1. Singular curves of the map 



curve for unequal-amplitude initial conditions 
of \^ with singular curves are marked by letters 'A, B, C, D, E 



\2.14i for negative e, as well as the initial-value 
2.21\) with e = —0.001. Five of the intersections 



can get, in addition to the vertical axis 70, are only two sequences 
{7(-l), 7(-l,-l), 7(-l,-l,-l), •••} 

and 

{7(1), 7(1,1), 7(1,1,1), ...}, 

which are shown in Fig. 16.11 Obviously these two sequences can not develop fractals. 
If we draw the initial-value curves (|4.ip and (14. 3p in this same plane, these curves 
will intersect with only a finite number of singular curves. An example is shown in 
Fig. 16.11 where the unequal-amplitude initial conditions (|2.2ip and e — —0.001 are 
taken. Thus the corresponding exit-velocity graph in the ipo space will have a limited 
number of singularity peaks, which is precisely what we observed in Fig. 11 (5)- (6) of 
Rcf. [25^. 

Now we quantitatively compare the map's predictions with the ODE results for 
negative values of e (comparison with the PDE results is expected to be similar, see 
Figs. I2.ll and l2.2p . For this purpose, we take e — —0.001 and the unequal-amplitude 
initial conditions (P?^ for the ODEs In the PDE (1^ with cubic-quintic 

nonlinearity (|2.15p . these ODE initial conditions correspond to S ~ —0.0003 and 
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unequal-amplitude initial conditions (j2.17p [25j . From the direct simulations of these 
ODEs, the graph of exit-velocity |C|oo versus the initial-phase difference ■00 is shown 
in Fig. I6.2f a). A finite number of singularity peaks can be seen, and this graph does 
not have a fractal structure. To obtain the map's predictions, we use the initial-value 
curve (|4.3p for this case as shown in Fig. 16.11 For each point on this initial-value 
curve, we iterate the map (|2.13p - (|2.14p to infinity to obtain the exit velocity from Eq. 
p.Sp . The results of the map's predictions are shown in Fig. 16.2( b). As can be seen, 
the map's prediction agrees with the ODE results very well. To better understand 
these graphs, we label five representative singularity peaks by letters 'A, B, C, D, E' 
in Fig. I6.2r b). Their corresponding points on the initial- value curve A*^ of Fig. 16. II are 
also labeled by the same letters. This connection makes it very easy to understand 
the exit-velocity graphs of Fig. 16.21 In particular, since has 10 intersections with 
singular curves in Fig. 16.11 this explains why the exit-velocity graphs in Fig. 16.21 
have 10 singularity peaks as well. One can further predict the dynamics of the ODE 
solution at each ipQ value in the exit-velocity graph of Fig. I6.2r a) in the same way 
as we did in the previous section. For instance, if ipo is on a hill of binary sequence 
a — (—1, —1, —1) in the exit-velocity fractal of Fig. 16.2( a). then the ({t) solution will 
oscillate three times before the collision time. The —ip solution will exhibit spikes at 
the three bounce points whose sign sequence is a. After the collision time, the (^(t) 
solution will go to —oo, while the tp{t) solution will approach a constant whose sign 
depends on which side of the hill the i/'o value is on. With this knowledge on the 
ODE solution, we can then predict the dynamics of the PDE solution, again in the 
same way as what we did in the previous section. If the two waves initially have equal 
amplitude, where the initial conditions for the ODEs are (|2.19p . then the initial- value 
curve A^ in (j4.ip will only intersect one singular curve 70 (i.e. the vertical axis). In 
this case, the exit-velocity graph of the ODE will have a unique singularity point at 
V'o = and be smooth elsewhere. The corresponding dynamics of ODE and PDE 
solutions can be predicted in the same way as above. With these results, an intimate 
knowledge is then obtained for the wave interactions in the e < case. 

7. Conclusion and discussion. In this paper, we analyzed the separatrix map 
which governs weak interactions of solitary waves in the generalized NLS equations. 
We showed that when e > 0, this map exhibits a fractal structure which we delineated 
by tracking its singular curves. Then through this map's fractal as well as connections 
between the map and the ODEs/PDEs, we reached a very deep understanding on the 
fractal structures as well as their interaction dynamics in the ODEs and PDEs. In 
addition, we analytically determined how the ODE and PDE fractals change as the 
parameter e varies, and showed that these predictions agree well with the numerical 
results. Furthermore, we proved a claim made earlier in |25j that fractal structures 
in the ODEs/PDEs for e > bifurcate from the singularity points in the integrable 
ODEs. When e < 0, we showed that the separatrix map does not possess a fractal 
structure, hence fractal scatterings can not occur in the corresponding ODEs and 
PDEs. 

Gathering all our results, we can now give a precise criterion for the existence of 
fractal scatterings in weak interactions of solitary waves in the generalized nonlinear 
Schrodinger equations. When \e\ ^ 1, fractal scatterings will occur if and only if the 
following two conditions are met: 

1. e > 0; 

2. The map's initial-condition curve A^ [see (|4.ip and (|4.3p for instance], under 
the restriction of iJ < 0, intersects the map's fractal region (see Figs. 13.11 14.11 
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Fig. 6.2. Comparison of the exit-velocity graphs from (a) ODE predictions, and (b) map 
predictions for e = —0.001. The initial conditions of the ODEs are If 2. 21}) for unequal- amplitude 
initial waves. Letters 'A, B, are locations of some singularity peaks in (b), which correspond to 
the intersections of \^ with singular curves in Fig. Ig.il (marked by the same letters). 



and[01). 

When this criterion is met (i.e. fractal scatterings occur), we can accurately predict 
the fractal structure as well as the interaction dynamics in the exit-velocity graph by 
simply drawing the initial-condition curve on the map's |rnoo| fractal of Fig. 13.11 
Thus, by now we have provided a simple answer to a complicated fractal-scattering 
problem in weak wave interactions. 

The only discrepancy between the PDE's exit-velocity graph and our ODE/map 
predictions is at the middle part of each hill, where the ODE's and map's graphs ex- 
hibit peaks of infinite height, but the PDE's graph dips down instead (see Figs. I2.1l and 
12. 2p . In those parameter regions, the two waves come together and interact strongly, 
which makes our reduced ODE model invalid. We have performed preliminary nu- 
merical investigations of the PDEs at the middle part of each hill in the exit- velocity 
graph, and found that the dips in those regions are not smooth. Inside each dip, we 
found finer structures which arc also fractal-like! These fractal-like structures inside 
each dip are apparently the product of strong wave interactions, thus should be re- 
lated to fractal scatterings in solitary wave collisions as reported in [TTJ [121 HSl [H] . 
Detailed investigations of finer structures inside dips of the PDE's exit- velocity graph 
lie outside the scope of this paper, and will be left for future studies. 
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